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Abstract 

We study the thermal properties of 54 Fe with the Brown-Richter interaction 
in the complete lpOf model space. Monte Carlo calculations show a peak 
in the heat capacity and rapid increases in both the moment of inertia and 
Ml strength near a temperature of 1.1 MeV that are associated with the 
vanishing of proton-proton and neutron-neutron monopole pair correlations; 
neutron-proton correlations persist to higher temperatures. Our results are 
consistent with a Fermi gas level density whose back-shift vanishes with in- 
creasing temperature. 

PACS numbers: 21.60.Cs, 21.60.Ka, 27.40,+z, 21.10.Ma 
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The nuclear level density increases rapidly at excitation energies above several MeV 
and it becomes difficult to resolve or calculate individual states. In this regime, it is more 
appropriate to employ a statistical description where observables are averaged over the many 
levels at a given energy. The concept of an equilibrated compound nucleus is among the most 
fundamental of nuclear reaction theories KJ, an d plays a central role in our understanding 
of processes induced by probes ranging from photons to heavy ions. 

While the proper description of a compound nucleus is in terms of a microcanonical 
(fixed-energy) ensemble, it is often more convenient to consider a canonical ensemble whose 
temperature is chosen to reproduce the average excitation energy. In the past decade, there 
has been renewed experimental and theoretical [[J effort to explore the properties of heavy 
nuclei at finite temperature and high spin. The properties of hot nuclei are also important 
in various astrophysical scenarios, particularly in the late stage of a supernova collapse and 
explosion ||. 

Most theoretical approaches to hot nuclei devolve to a mean-field description based on 
an average configuration ||. The realization that thermal and quantal fluctuations about 
the average can be important has prompted more sophisticated approximations || , although 
even these have clear limitations. In principle, the nuclear shell model (which provides a 
complete spectrum and wavefunctions) offers a fully microscopic approach to the problem. 
However, conventional finite-temperature shell model calculations within a complete major 
shell are limited to light nuclei ( 20 Ne and 24 Mg) in the sd shell 0. 

In this Letter, we exploit recently developed Monte Carlo techniques to calculate the 
thermal properties of 54 Fe in a complete Ohu model space with a realistic interaction. The 
methods we use describe the nucleus by a canonical ensemble at temperature T = and 
employ a Hubbard-Stratonovich linearization of the imaginary-time many-body propaga- 
tor, e~@ H , to express observables as path integrals of one-body propagators in fluctuating 
auxiliary fields ||. Since Monte Carlo techniques avoid an explicit enumeration of the many- 
body states, they can be used in model spaces far larger than those accessible to conventional 
methods. The Monte Carlo results are in principle exact and are in practice subject only to 
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controllable sampling and discretization errors. The nucleus we have chosen for this initial 
study ( 54 Fe) is among the most abundant in the presupernova core of a massive star, so that 
its thermal properties are of considerable astrophysical import. Further, the ground states 
of nuclei in the mid-p/-shell are dominated by nucleon-nucleon correlations (e.g., pairing) 
whose evolution with increasing temperature is of particular interest. 

To circumvent the "sign problem" encountered in the Monte Carlo shell model calcu- 
lations with realistic interactions, Alhassid et al. || suggested an extrapolation procedure 
from a family of Hamiltonians that are free of the sign problem to the physical Hamiltonian. 
One defines a set of Hamiltonians H g = Hq + qHb such that H g= i = H is the physical 
Hamiltonian and Hq,b are the "good" and "bad" parts of the Hamiltonian, respectively. 
For g < 0, H g is free of the sign problem and calculated observables are extrapolated to 
g = 1. For ground state properties, this procedure was validated by comparison to direct 
diagonalization results in the sd and lower pf shells. However, it is impractical at inter- 
mediate temperatures due to the overly strong pairing interaction in H g for g < 0, which 
suppresses the population of excited states. This problem can be corrected by scaling Hq 
as (l — Hq, together with a (^-dependent compression of the single-particle spectrum; 
the value of x is chosen to make the ^-extrapolation as smooth as possible. Note that, as 
before, the original Hamiltonian is recovered for g = 1. 

Our calculations include the complete set of 1^3/2,1/20/7/2,5/2 states interacting through 
the realistic Brown- Richter Hamiltonian |10|| . Some 5 x 10 9 configurations of the 8 valence 



neutrons and 6 valence protons moving in these 20 orbitals are involved in the canonical 
ensemble. The results presented below have been obtained in MC shell model studies with 
a time step of A/3 = 1/32 MeV -1 using 5000-9000 independent Monte Carlo samples at 
seven values of the coupling constant g spaced between —1 and and the value \ — 4. 
A linear extrapolation to the physical case (g = 1) is justified by the quality-of-fit for 
most of the observables discussed below, although the quadrupole moments warranted a 
quadratic extrapolation. We have tested our procedure for the /p-shell nucleus 44 Ti against 
a calculation performed with the direct diagonalization code CRUNCHER [|ll|] and were able 



to reproduce the excitation energy as a function of temperature for the temperature interval 
relevant for this paper ||12j| . However, we note that an exact reproduction of the CRUNCHER 
energies is obtained only after a A/3 — > extrapolation, which lowers the absolute energies 
slightly compared to the calculation with the finite value A/3 = 1/32 MeV -1 . For 54 Fe we 
have checked at several temperatures that our qualitative results are not changed by the 
A/3 extrapolation. We have also checked that using x = 3 in the ^-extrapolation does not 
change our results. 

The calculated temperature dependence of various observables is shown in Fig. |l|. In 
accord with general thermodynamic principles, the internal energy U steadily increases with 
increasing temperature [pl|. It shows an inflection point around T pa 1.1 MeV, leading to 



a peak in the heat capacity, C = dU / dT, whose physical origin we will discuss below. The 
decrease in C for T > 1.4 MeV is due to our finite model space (the Schottky effect [0); 
we estimate that limitation of the model space to only the p/-shell renders our calculations 
of 54 Fe quantitatively unreliable for temperatures above this value (internal energies U ^ 
15 MeV). The same behavior is apparent in the level density parameter, a = C/2T. The 
empirical value for a is A/8 MeV = 6.8 MeV" 1 which is in good agreement with our results 
for T pa 1.1-1.5 MeV. 

We also show in Fig. [I] the expectation values of the squares of the J = proton- 
proton and neutron-neutron pairing fields, (A* A). Although the pair wave function we 
have used (the BCS form, in which all time-reversed pairs have equal amplitudes) is not 
necessarily optimal, these observables are a rough measure of the number of J = pairs 
in the nucleus. At low temperatures, the pairing fields are significantly larger than those 
calculated for a non-interacting Fermi gas, indicating a strong coherence in the ground state. 
With increasing temperature, the pairing fields decrease; both approach the Fermi gas values 
for T pa 1.5 MeV and follow it closely for even higher temperatures. Associated with the 
breaking of pairs is a dramatic increase in the moment of inertia, J, for T = 1.0-1.5 MeV; 
this is analogous to the rapid increase in magnetic susceptibility in a superconductor. At 
temperatures above 1.5 MeV, I is in agreement with the rigid rotor value, 10.77i 2 /MeV; at 



even higher temperatures it decreases linearly due to our finite model space. 

In Fig. |2|, we show various static observables. The Ml strength unquenches rapidly with 
heating near the transition temperature. However, for T = 1.3-2 MeV B(Ml) remains sig- 
nificantly lower than the single-particle estimate (41 fj,%) , suggesting a persistent quenching 
at temperatures above the like-nucleon depairing. This finding is supported by the near- 
constancy of the Gamow- Teller (5 + strength, B(GT + ), for temperatures up to 2 MeV. As the 
results of Ref. fl5j demonstrate that neutron-proton correlations are responsible for much 



of the GT quenching in iron nuclei at zero temperature, we interpret the present results as 
evidence that isovector proton-neutron correlations persist to higher temperatures. We have 
verified that, in our restricted model space, both the GT + and Ml strengths unquench at 
temperatures above 2 MeV and that, in the high-temperature limit, they both approach the 
appropriate Fermi gas values. We note that it is often assumed in astrophysical calculations 
that the GT strength is independent of temperature ||; our calculations demonstrate that 
this is true for the relevant temperature regime (T < 2 MeV). We also note that a detailed 
examination of the occupation numbers of the various orbitals show no unusual variation as 
the pairing vanishes. 

The isoscalar mass quadupole moment, (Q 2 ), increases at temperatures near the phase 
transition (Fig. ||, upper right), while the isovector moment, (Ql), decreases as the nucleus is 
heated, showing a minimum near T = 1.1 MeV (Fig. ^ lower right). The behaviors of (Q 2 ) 
and {Q 2 ) imply that {Q p ■ Q n ) = {{Q 2 ) — {Q 2 ))/^ increases near the phase transition (Q Pjn 
are the proton and neutron quadrupole moments). Noting the relation between {Q p -Q n ) and 
the orbital part of the Ml strength ||16|| , we interpret the partial unquenching of the B (Ml) 



near T = 1.1 MeV as related to the orbital part, while the spin part, dominated by the same 
operator as the GT strength, remains significantly quenched to higher temperatures. 

We have compared our results for U with two simple models. As in Ref. [[IT]], we define 
the partition function of the nucleus as 

roc 

Z = V(2J t + l)e~ pEi + / p{E)e- f3E dE , (1) 
i Je 



where the sum runs over the experimentally known nuclear levels and the continuum state 
density p(E) has been approximated by the backshifted Fermi gas model |I5| , where the 
backshift P accounts for the energy to break a pair; we adopted a level density parameter 
a = 7.2 MeV -1 and chose Eq = 4 MeV to smoothly match the two terms in Eq. (1). The 
first model had the conventional T-independent backshift Pq = 1.45 MeV [|18|| , while the 
second model simulated the vanishing of the pairing by a temperature-dependent P: 

P(T)=P (l + exp{^°}) 1 . (2) 

Our choice of the parameters To = 1.05 MeV, a = 0.25 MeV was motivated from the two 
upper right panels in Fig. [I]. As seen from the solid curves in the left panels of Fig. [I], 
the P{T) model is in better agreement with the Monte Carlo shell model results than is 
the constant-P model; in particular, there are clear maxima in both the heat capacity and 
the level density parameter related to the pairing phase transition. The dashed curves in 
these panels indicate that the assumption of a constant level density parameter is not quite 
appropriate at temperatures below the phase transition. 

In conclusion, we have demonstrated that shell model Monte Carlo methods are well- 
suited to studying the finite-temperature properties of nuclei using realistic effective two- 
body interactions. Our calculations of 54 Fe in the complete p/-shell show clear signatures of 
a pairing phase transformation at a temperature of 1.1 MeV, but persistent quenching of the 
Gamow- Teller f3 + strength at higher temperatures. Results at temperatures above 1.5 MeV 
will become reliable only when two or more major shells are included in the calculations, 
an elaboration that is computationally quite feasible. The extension of these calculations 
to other interactions, heavier nuclei, and other observables should allow a more thorough 
understanding of nuclear properties at high excitation energies than is now possible by other 
methods. 
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FIGURES 

FIG. 1. Temperature dependence of various observables in 54 Fe. Monte Carlo points with sta- 
tistical errors are shown at each temperature T. In the left-hand column, the internal energy, U, is 
calculated as {H)—Eq, where H is the many-body Hamiltonian and Eq the ground state energy. The 
heat capacity C is calculated by a finite-difference approximation to dU/dT, after U(T) has been 
subjected to a three-point smoothing, and the level density parameter is a = C/2T. The dashed 
and solid curves in these panels correspond to the constant- and temperature-dependent-backshift 
Fermi gas models, as described in the text. To eliminate the systematic error associated with the 
determination of Eq, we have chosen this parameter so that the Monte Carlo and Fermi gas results 
for U are equal at T = 0.66 MeV. In the right-hand column, we show the expectation values of the 
squares of the proton and neutron pairing fields, where = J2Pj m Pjm ( an d similarly for neu- 
trons) and the sum is over all orbitals with m > 0. For comparison, the pairing fields calculated in 
an uncorrelated Fermi gas are shown by the solid curve. The moment of inertia is obtained from 
the expectation values of the square of the total angular momentum by / = (3(J 2 )/3. 

FIG. 2. The upper left panel shows the total magnetic dipole strength, B(M1) in units of 
nuclear magnetons; it is calculated using free-nucleon g- factors. The lower left panel shows the GT 
/3+ strength. The low-temperature value of 4.7 ± 0.2 is somewhat larger than the 4.3 ± 0.2 given 
in Refs. [9,15], since the latter calculations were done with A/3 = 1/16 MeV -1 . Calculations at 
A/3 = 1/64 MeV -1 indicate that the present low-temperature value is converged. The upper and 
lower right panels show the isoscalar (Q = Q p + Q n ) and isovector (Q v = Q p — Q n ) quadrupole 
strengths, where the quadrupole operators are r 2 Y2 and results are given in terms of the oscillator 
length, b = 1.96 fm. 
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